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16.  Akstrwt 


A 


This  report  describes  an  efficient,  autocorrelation  based  algorithm  for  estimating  low  altitude  radial 
winds  u^ing  signals  from  the  two  receiving  beams  of  an  airport  surveillance  radar  (ASR).  The  approach 
seeks  to  achieve  the  accuracy  demonstrated  previously  for  spectral  domain  dual  beam  velocity  estimators 
with  significantly  reduced  computational  requirements.  Fundamental  to  the  technique  is  the  assumption 
that  the  power  spectrum  measured  with  an  airport  surveillance  radar's  broad  elevation  beam  can  be  fitted 
by  a  two  component  Gaussian  model.  The  parameters  of  this  model  are  estimated  using  measured  low- 
order  autocorrelation  lags  from  the  low  and  high  beam  received  signals.  The  desired  near  surface  radial 
velocity  estimate  is  obtained  directly  as  one  of  these  parameters  —  the  center  frequency  of  ihef^low  ^ 

allilud^lGau.ssian  spectrum  component.  A.  v'  // 

Simulated  data  and  field  measurements  from  l.incoln  l.ab64»U»F;J’s  experimental  ASR-8  in  Huntsville,  A. 
•Abibama' were  used  to  evaluate  the  accuracy  of  the  autocorrelation  based  velocity  estimates.  Monto  (iarlo 
simulations  indicate  that  biases  relative  to  the  near  surface  outflow  velocity  in  a  microhursi  would  be  less 
than  2.5  m/s  unless  the  microburst  were  distant  (range  >  12  km)  or  very  shallow  (depth  of  maximum  wind 
speed  layer  <  50  m).  Estimate  standard  deviations  averaged  0.5  m/s  after  the  spatial  filtering  employed  in 
our  processing  sequence.  The  algorithm's  velocity  estimate  accuracy  was  sufficient  to  allow  for  automatic 
detection  of  measured  microbursts  during  1988  with  a  detection  probability  exceeding  0.9  and  a  false 
alarm  probability  le.ss  than  0.05.. Our  analysis  indicates  that  the  dual-beam  autocorrelation  based  velocity 
estimator  should  support  ASR  wind  shear  detection  at  approximately  the  same  level  of  confidence  as  the 
low-high  beam  spectral  differencing  algorithm  evaluated  by  Weber  ami  Noyes  (1088).  ^  _ 
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Dual-Beam  Autocorrelation  Based  Wind  Estimates  from 
Airport  Surveillance  Radar  Signals 


ABSTRACT 

This  report  describes  an  efficient,  autocorrelation  bfised  algo¬ 
rithm  for  estimating  low  altitude  radial  winds  using  signals  from  the 
two  receiving  beams  of  an  airport  surveillance  radar  (ASR).  The 
approach  seeks  to  achieve  the  accuracy  demonstrated  previously  for 
spectral  domain  dual  beam  velocity  estimators  with  significantly 
reduced  computational  requirements.  Fundamental  to  the  technique 
is  the  assumption  that  the  power  spectrum  measured  with  an  air¬ 
port  surveillance  radar’s  broad  elevation  beam  can  be  fitted  by  a 
two  component  Gaussian  model.  The  parameters  of  this  model  are 
estimated  using  measured  low-order  autocorrelation  lags  from  the 
low  and  high  beam  received  signals.  I’he  desired  near  surface  radial 
velocity  estimate  is  obtained  directly  as  one  of  these  parameters  — 
the  center  frequency  of  the  "low  altitude"  Gaussian  spectrum  com¬ 
ponent. 

Simulated  data  and  field  measurements  from  Lincoln 
Laboratory’s  experimental  ASR-S  in  Huntsville*,  Alabama  were  used 
to  evaluate  the  accuracy  of  the  autocorrelation  based  velocity  esti¬ 
mates.  Monte  Carlo  simulations  indicate  that  biases  relative  to  the 
near  surface  outflow  velocity  in  a  microburst  would  be  less  than  *2.5 
m/s  unless  the  microburst  were  distant  (range  >  12  km)  or  very 
shallow  (depth  of  maximum  wind  speed  layer  <  50  m).  Estimate 
standard  deviations  averaged  0.5  m/s  after  the  spatial  filtering 
employed  in  our  proce*ssing  sequence.  The  algorithm’s  velocity  esti¬ 
mate  accuracy  was  sufficient  to  allow  for  automatic  detection  of 
measured  microbursts  during  1988  with  a  detection  probability 
exceeding  0.9  and  a  false  alarm  probability  less  than  t).05.  Our 
analysis  indicates  that  the  dual-beam  autocorrelation  based  velocity 
estimator  should  siqiport  ASR  wind  shear  detection  at  approxi¬ 
mately  the  same  level  of  confidence  as  the  low-high  lieam  spectral 
differencing  algorithm  evaluated  by  Weber  and  Noyes  (1988). 
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Images  of  the  radial  velocity  field  in  a  inicrol)urs(  iiroduciiig 
thunderstorm  near  Huntsville,  Alabama  on  Ih  August  1988. 
The  upper  left  panel  shows  measurements  from  the  MiT  C- 
band  beam  weather  radar  scanning  in  PPl  mode  at  O.fi" 
elevation  angle.  Upper  right  panel  is  lield  estimated  from 
ASR  signals  using  the  spectral  (liU'erencing  method.  Lower 
left  panel  is  corresponding  field  from  autocorrelation  based 
algorithm.  Lower  right  panel  nsre  autocorrelation  based 
algorithm  and  high  and  low  beam  signals  collected  on 
different  antenna  .scans. 

Images  of  the  radial  velocity  lield  in  a  inicroburst  lu'oducing 
thunderstorm  near  lluntsvillc,  .Alabama  on  2t  .June  1988. 
The  upper  panel  shows  mea.surements  from  the  .MIT  C'-band 
beam  weather  radar  scanning  in  PIM  mode  at  0.6"  elevation 
angle.  Lower  right  panel  is  field  estimated  frcmi  A.SR  signals 
using  the  spectral  differencing  method.  Lower  left  panel  is 
corresponding  lield  from  autocorrelation  based  algorithm. 
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thunderstorm  near  Huntsville.  Alabama  on  21  .June  1988. 
d'he  upper  panel  shows  measurements  from  the  MPf  C-band 
b('am  weather  radar  scanning  in  I*P1  mode-  at  0.6"  elevation 
angle.  Lower  right  panel  is  lield  estimated  fi-om  .\SR  signals 
using  the  spectral  differencing  method.  Lower  left  imnel  is 
corresponding  field  from  autocorrelation  based  algorithm. 

linages  of  j-adial  velocity  field  i))  a.  ))jicrobm‘st  pj-oducing 
thunderstorm  near  Huntsville.  .Alabama  on  2S  .June  1988. 
The  upper  panel  shows  measurements  from  the  MIT  C-band 
lieam  weather  radar  scanning  in  I’I’l  nu.xle  ;it  0.6"  elevation 
angle.  Lower  right  i)anel  is  field  estimated  from  .ASR  signals 
using  the  spectral  dillereucing  method.  Lower  left  panel  is 
corresponding  field  from  autocorrelation  based  algorithm. 

Simulated  underlying  .ASR  |)()wer  sj)ecti-a  foi-  microl)urst  wind 
field  model  of  Figure  I1I-2  (and  its  miiror  image).  I^lot  for¬ 
mat  is  as  in  Figures  If-I  and  il-2. 

Simulated  three-degree-of-freedom  .AtsR  power  s|)ectra  esti¬ 
mates.  I Jnderlying  spectra  are  as  in  l''igure  C-l. 

Spectral  differencing  and  autocorrelation  based  velocity  esti¬ 
mate  bias  and  standard  deviation  versus  "true"  outflow  velo¬ 
city.  Tlie  simulations  used  hen'  r('produee  the  alternating 
I’RF  waveform  of  an  .ASR-9  and  the  use  of  a  four-eoc'flicient 
interpolation  filter.  I’lot  formal  is  as  in  I'iguia'  111-1. 
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Dual-Beam  Autocorrelation  Based  Wind  Estimates  from 
Airport  Surveillance  Radar  Signals 


I.  INTRODUCTION 

This  report  evaluates  an  efficient,  autocorrelation  based  algoritlmi  for  estimat¬ 
ing  low  altitude  radial  winds  using  signals  from  the  two  receiving  beams  of  an  air¬ 
port  surveillance  radar  (ASR).  The  described  approach  seeks  to  achieve  the  accu¬ 
racy  demonstrated  previously  for  spectral  domain  dual  beam  velocity  estimators 
(Anderson, 1989;  Atlas, 1987, 1989;  Weber  and  Noyes, 1988)  with  significantly 
reduced  computational  requirements. 

The  feasibility  of  a  data  processing  augmentation  to  existing  ASRs  that  would 
allow  for  detection  of  low  altitude  wind  shear  has  been  under  investigation  by 
Lincoln  Laboratory  and  cooperating  universities  since  1984.  This  capability  would 
require; 

(a)  a  signal  processing  module  that  would  suppress  ground  clutter  and  estimate 
the  radial  velocity  of  near-surjace  precipitation  wind  tracers; 

(b)  an  algorithm  to  automatically  recognize  significant  wind  shear  in  the  result¬ 
ing  velocity  field. 

.\nderson  (1987),  Weber  and  Moser  (1987)  and  Weber  (1987)  considered  the 
impact  of  ground  clutter  on  ASR  wind  measurements.  Simulations  and  analysis 
of  real  clutter  data  from  airport  surveillance  radars  indicated  that  wind  measure¬ 
ments  were  feasible  even  at  short  range  provided  that  the  reflectivity  factor  of  the 
precipitation  tracers  was  approximately  20  dBz  or  greater.  Even  when  this  condi¬ 
tion  is  met  however,  accurate  low-altitude  radial  wind  estimates  may  not  be 
obtained  with  "conventional"  mean  Doppler  estimators  when  the  radial  wind  com¬ 
ponent  varies  rapidly  with  altitude.  In  this  situation,  an  ASR's  fan-shaped  eleva¬ 
tion  beam  intercepts  scatterers  moving  at  both  the  near-surface  velocity  and  the 
velocity  of  wind  aloft;  the  result  is  a  broad,  possibly  multi-modal  velocity  spec¬ 
trum  whose  power-weighted  mean  differs  markedly  fnnn  the  near  surface  radial 
wind  velocity.  Since  both  microbursts  and  gust  fronts  exhibit  strotig  vertical 
shear  in  the  horizontal  winds  near  the  grouiul,  this  beam  resolution  issue  is  of 
obvious  importance. 

Weber  and  Noyes  (1988)  used  data  from  an  experimental  ASR  operattsl  during 
periods  of  nearby  thiinderstcM’in  activity  to  evaluate  three  methods  for  ('stimating 
low  altitude  winds  from  ASR  signals; 

(i)  high-pass  filtering  to  exploit  the  fact  that  microbiirst  oiitilow  winds  ai-e  oft en 
higher  in  al.'scdiite  magnitude  than  winds  aloft; 

(ii)  comparison  of  the  strength  of  divergence  regions  d('t('cl(Ml  in  velocity  fields 
from  the  high  and  low  |■eceiving  beams  of  the  ASR  to  "coi'rect"  the  measured 
velocity  shear  values; 

(iii)  resolution  cell  by  cell  com|)aris.on  of  the  |)ow(a'  s|>ectra  of  the  low  and  high 
receiving  la'ams  to  det (’|•min('  the  velocity  domain  associated  with  lU'ar  sur¬ 
face  scatteix'rs. 
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Our  assessment  included  an  end-toend  data  processing  evaluation  whereby  wind 
fields  estimated  from  the  ASR  signals  were  passed  into  a  slightly  modified  version 
of  the  Terminal  Doppler  Weather  Radar's  (TDWIi)  surface  outtlow  detection  algo¬ 
rithm  (Merritt,  1987).  The  resulting  microbuist  alarms  were  scored  on  a  scan-by¬ 
scan  basis  against  "truth"  as  determined  from  manual  examination  of  wind  fields 
measured  with  a  colocated  pencil-beam  Doppler  weather  radar.  Overall,  the  third 
method  above  —  transformation  of  high  and  low  beam  signals  into  tlie  frequency 
domain  followed  by  spectral  difl’erencing  —  produced  better  detection  performance 
(higher  probability  of  detection  and  tower  false  alarm  j)robability)  and  more  accu¬ 
rate  velocity  shear  estimates  than  the  other  two  methods  evaluated. 

Use  of  this  dual-beam  signal  processing  approach,  however,  raises  several 
implementation  issues.  ASRs  normally  transmit  circularly  polarized  (C’i*)  signals 
during  heavy  rain  in  order  to  reduce  precipitation  clutter  in  the  aircraft  detection 
processor.  In  order  to  avoid  a  resulting  l.'i  to  20  dR  loss  in  received  power  from 
weather  echoes,  any  weather  processor  on  an  airport  surveillance  radar  should 
receive  its  input  from  the  opposite-sense  |)olarizcd  antenna  i)ort.  ASR-8's  and 
ASR-9’s  have  only  one  path  through  the  rotary  joint  for  opposite  sense  polarized 
signals:  thus  during  operation  with  circular  polarization,  weather  data  from  both 
high  and  low  beams  could  not  be  accessed  simultaneously.  Rower  spectra  for  the 
two  beams  would  have  to  be  calculated  on  alternate  antenna  scans,  requiring 
memory  storage  for  all  data  from  one  of  the  scans.  Assuming  range  coverage  to 
20  km  and  two  byte  integer  representation  of  the  in-pliase  and  (piadrature  signals, 
this  would  require  3  Mbytes  of  dedicated  physical  imunory. 

Another  issue  is  the  conqnjtational  re(|uirem('nl  of  ('stimating  power  sf)ectra  for 
the  two  beams  in  each  resolution  cell  and  calculating  a  velocity  based  on  the 
difference  of  the  two  spectra.  Assuming  that  the  spectra  were  estimated  using 
Fast  Fourier  Transforms  (FFT),  about  lOOO  Hoating  point  operations  per  resolu¬ 
tion  cell  would  be  rccjuired.  For  the  same  20  km  range  coverage  considered  prevn- 
ously,  this  translates  to  17  million  Hoating  point  operations  per  secoiul.  While 
l)oth  the  memory  storage  and  processing  speed  re(|uirements  could  be  met.  they 
w'ould  certainly  drive  the  cost  of  the  signal  lu-ocessing  computer  higher. 

A  third  issue  arises  from  the  large  variance  and  coarse  (luantization  of  power 
spectrum  estimates  obtained  from  a  rapidly  scanning  aiiqiort  surveillance  radar. 
As  implemented  by  Weber  and  Noyes  (1988).  the  low  and  high  Ixavin  power  spec¬ 
trum  estimates  were  calculated  from  31  successive'  jeulses  (2.6®  in  azimuth)  and 
incoherently  averaged  over  three  successive  range  gates  (3(50  m).  Velocity  n'solu- 
tion  wa.s  about  2  m/s  and  the  90*^7  conlidence  interval  of  the  spectral  estimates 
was  7  (in.  Monte  C'arlo  simulations  indicate  that  the  stamlard  deviation  of  mean 
velocity  estimates  using  our  implementation  is  a|)proximately  2  m/s.  'Phus  addi¬ 
tional  spatial  filtering  had  to  be  applied  to  tin'  velocity  field  in  order  to  achieve 
acceptable  yx'rformance  from  the  inicroburst  deti'ction  algorithm. 

In  this  re|)ort.  we  evaluate  an  alt('rnativ('  dual-beam  velocity  estimation  tech- 
ni<|ne  that  significantly  reduces  both  memory  storage  and  computational  reejnire- 
nients.  'Phe  basic  strategy  is  to  estimate  the  i)arani('ters  of  an  assumed  biinodal 
(laiissian  power  spectra,  based  on  low-ord('r  autocorrelation  lags  from  the  hnv  and 
high  beam  received  signals,  d'he  desin'd  low-altitiuh'  vc'locity  estimate  is  then 
obtained  directly  as  the  center  freciuency  of  one  of  the  two  components  of  this 
s|)ectral  model.  Ih'cause  only  the  zero  and  one  sami)le  delay  autocorrelation  lags 
for  each  beam  are  \ised.  |)hysical  memory  storage  re(iuirements  would  be  reduced 
to  ().')  Mbyt('.  C'om|nit at ional  re(|uireinents  art'  alxnit  2  million  Hoating  jioint 


operations  per  second.  Estimate  standard  deviations  of  2  m/s  can  be  achieved 
without  the  need  for  range  smoothing  as  required  for  the  spectral  differencing 
approach. 

Section  II  describes  the  double  Gaussian  spectral  model  and  algorithms  for 
estimating  its  parameters  from  autocorrelation  function  measurements.  This  dis¬ 
cussion  draws  heavily  on  measured  microburst  power  spectra  from  our  experimen¬ 
tal  ASR.  In  Section  III,  we  evaluate  the  performance  of  the  algorithm,  using 
simulated  and  measured  ASR  signals  from  microbursts.  Automatic  microburst 
detection  performance  using  this  method  and  the  spectral  differencing  approach 
evaluated  by  Weber  and  Noyes  (1988)  is  compared  statistically  for  our  1988  field 
measurements.  Over  the  evaluated  data  set,  the  two  methods  provide  similar  per¬ 
formance.  Section  IV  summarizes  our  findings  and  describes  ongoing  work. 


n.  VELOCITY  ESTIMATION  BASED  ON  A  DUAL  GAUSSIAN  SPEC¬ 
TRUM  SHAPE  MODEL 


A.  Velocity  Spectra  In  Microbiirst  Cores 

Figure  II- 1  shows  examples  of  velocity  s[)ectra  measured  with  our  experimental 
airport  surveillance  radar  in  the  cores  of  eight  "wet"  microbursts  near  Huntsville, 
Alabama.  The  spectra  were  estimated  from  hTTs  of  34-sample  Hamming  win¬ 
dowed  data  sequences,  incoherently  averaged  over  three  contiguous  range  gates. 
An  adaptively-selected  clutter  filter  (Weber,  1987)  was  applied  to  the  signals  prior 
to  F'ourier  transformation;  in  most  of  the  examples  shown,  however,  the 
reflectivity  factor  in  the  microburst  core  was  sufficiently  high  that  no  filtering  was 
required.  Although  our  experimental  radar  transmitted  pulses  at  a  uniform 
repetition  frequency  (PRF),  the  eight/ten-pulse  alternating  PRF'  sequence  to  be 
used  by  ASR-9s  could  be  accommodated  by  preceding  the  clutter  filters  with  a 
shift-variant  interpolation  filter  as  described  by  Weber  (1987). 

The  plots  show  calculated  spectra  for  both  high  (dashed)  and  low  (solid)  receiv¬ 
ing  beams.  The  left  hand  panels  are  for  the  approaching  radial  velocity  cores  of 
the  microbui-sts  and  those  in  the  right  column  are  for  the  receding  cores.  The 
spectra  have  been  normalized  so  that  the  areas  under  the  curves  are  equal.  For 
reference,  radial  velocities  measured  by  the  pencil  beam  weather  radar  at  the  same 
range-azimuth  locations  and  times  are  imiicated  by  dashed  vertical  lines.  The 
pencil  beam  radar  was  scanned  at  O.ti"  elevation  angle  to  estimate  the  near  surface 
radial  wind  sp^-ed. 

As  discussed  in  Weber  and  Noyes  (1988)  the  broad  width  and  multimodal 
nature  of  these  power  spectra  result  in  significant  differences  between  the  "true" 
near  surface  radial  velocity  and  the  A.SR-based  measurement  when  a  conventional 
mean  Doppler  estimator  is  used.  The  discrepancy  results,  of  course  because  the 
applicability  of  the  mean  Doppler  estimate  rests  on  the  assumption  of  narrow- 
width,  roughly  symmetric  power  spectral  shaiie.  Of  the  displayed  power  spectra, 
only  those  in  the  approaching  core  on  I  August  exhibit  symmetry.  The  remaining 
sfiectra  are  strongly  skewed  and/or  hi-  or  multi-modal.  Some  of  this  c(.m])lexity 
arises  from  spectral  estimate  variance  (recall  that  our  jirocedure  results  in  a  7  dR 
9()^V  confidence  interval).  As  illustrated  in  Appendix  C  however,  the  dominant 
features  of  the  spectra  can  be  attributed  to  the  interaction  of  an  ASR’s  fan-shapc'd 
eh'vation  beam  with  the  strongly  sheared  radial  wind  field  in  an  a  microburst 
out  tlow. 

Note  that  the  low  beam  |)ovvpr  spectral  density  always  exceeds  that  of  the  high 
Ix'am  near  the  indicated  surface  w'ind  speed.  This  rellects  reduced  high  beam 
antenna  gain  near  the  surface',  and  is  the  basis  for  the  spectral  differencing  tech- 
ni(|ue  alluded  to  above. 
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[•'igiire  11-1.  Velocity  spectra  measured  with  ASR  in  approaching  (left)  and  reced¬ 
ing  (right)  radial  velocity  cores  of  example  microbursts.  The  ordinate  is  relative 
power  in  linear  units.  Spectra  are  plotted  for  both  low  (solid)  and  high  (dashed) 
l)eam  signals.  l)a.slied  vertical  lines  show  mean  velocities  measurefl  by  pencil 
beam  radar  at  0.7  degrees  elevation  at  same  locations  and  times. 


B.  Autocorrelation  Based  Dual  Gaussian  Fits  to  the  Measured  Spectra 
I.  Estimates  Using  R{t)  and  R{2t) 

Examination  of  the  above  measured  spectra  suggests  that  the  parameters  of  a 
double  Gaussian  power  spectrum  model  might  provide  a  better  estimate  of  the 
low  altitude  wind  field  than  the  conventional  mean  Doppler  method.  Physically, 
the  presence  of  two  separated  modes  in  the  power  spectra  implies  a  radial  velocity 
field  that  is  strongly  discontinuous  within  the  radar’s  beamwidth,  as  at  the  top  of 
a  shallow  microburst  outflow.  In  a  region  where  the  radial  wind  speed  varies 
linearly  with  height  or  is  constant,  the  double  Gaussian  model  can  also  adequately 
represent  the  power  spectrum  that  would  be  measured  by  an  ASR  since  the  two 
components  may  overlap  substantially,  producing  a  broad,  unimodal  spectrum.  \ 


We  assume  therefore  that  the  weather  echo  power  spectra  measured  by  the  fan 
shaped  elevation  beam  of  an  ASR  can  be  approximately  represented  as: 


^i{f  )  \/7r~  ^^1^1  «  I  "b 

V27r(7i 


> 


‘lar 


-ex  Pi 


2cr.? 


(1) 


where  i=l,2  Indexes  the  low  and  high  receiving  beams.  The  center  frequencies 
and  widths  of  the  two  spectral  components  are  assumed  to  be  identical  between 
the  two  beams  but  the  amplitudes  will  differ  owing  to  the  different  weightings 
from  the  antenna  patterns. 


Appendix  A  describes  a  method  of  estimating  the  eight  parameters  in  equations 
(1)  using  measurements  of  the  low  and  high  beam  autocorrelation  functions  at 
lags  T  and  2r.  The  solution  is  constructed  using  knowledge  of  the  beam  patterns 
so  that  the  parameters  a,  ,,  ctj  and  /,  represent  the  spectral  component  associated 
with  scatterers  at  low  elevation  angle.  Figure  II-2  plots  the  resulting  dual  Gaus¬ 
sian  power  spectra  for  the  microburst  cores  shown  previously.  Although  these  are 
not  optimum  dual  Gaussian  fits  to  the  measured  spectra  (for  example,  in  a  least 
squares  sense)  they  generally  correspond  well  to  the  data. 


Low  altitude  velocity  estimates  can  be  derived  from  the  calculated  spectral 
parameteis  by  reconstructing  the  spectra  and  implementing  a  low-high  beam 
differencing  algorithm  as  described  in  Weber  and  Noyes  (1988).  Table  Il-l  com¬ 
pares  the  resulting  velocity  shear  estimates  for  these  eight  microbui'sts  to  those 
measured  by  the  pencil  beam  "truth"  radar  and  to  those  calculated  using  the  spec¬ 
tral  differencing  technique  applied  directly  to  the  measured  power  spectra.  In 
seven  of  the  eight  cases  considered,  both  ASR-based  velocity  shear  estimates  are  in 
good  agreement  with  the  pencil  beam  measurements.  For  the  spectra  measured  on 
9-11-87,  dual  Gaussian  fitting  resulted  in  a  substantially  larger  wind  shear  esti¬ 
mate  than  given  by  the  pencil  beam  weather  radar  and  the  spectral  differencing 
technique;  this  produced  the  larger  overall  RMS  error  associated  with  (he  dual 
G  anssi  an  a j  >  p  rcjac  h . 


The  above  results  indicate  that  the  lidelity  of  the  dual  Gaussian  spectral  model 
to  the  measured  signal  characteristics  may  be  sufficient  to  generate  a  reliable  low 
altitude  velocity  estimate  from  ASR  data.  However,  the  described  solution  is  not 
computationally  efficient.  Calculation  of  the  first  and  second  autocorrelation  lags 
requires  l..'j  times  as  many  operations  as  would  be  needed  to  estimate  /f(())  and 
R{t).  More  imi)ortantly,  the  double  Gaussian  spectral  parameters  ar('  determined 
iteratively  (Appendix  A),  followed  by  numerical  integraticni  to  d('riv('  a  velocity 
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Figure  II-2.  Velocity  spectra  reconstructed  from  ASR  autocorrelation  lag  measure¬ 
ments  in  approaching  (left)  and  receding  (right)  radial  velocity  cores  of  micro¬ 
bursts.  The  spectral  model  of  equation  (1)  was  used  with  parameters  estimated 
from  /?(r)  and  R{2t)  as  described  in  Appendix  A.  The  microburst  examples  and 
plot  format  are  as  in  Figure  III-l. 


Power  (Linear  Units) 


Figure  II-2.  (continued) 


Table  II-l:  Microburst  Velocity  Shear  Estimates  Using  Spectral 
_ DifYerencintt  and  Rir).  R(2t]  Dual  Gaussian  Fit _ 1 _ 

Date/Time 

Pencil  Beam 

Radar  AV  (m/s) 

ASR  Specti-al 

Differencing 

(Measured) 

ASR  Spectral 
Dilferencing 
(Dual  Gaussian 
Fitl 

23 

21 

23 

5-21-87  14:14 

28 

28 

27 

6-21-87  20:45 

27 

29 

30 

6-21-87  20:38 

23 

23 

22 

9-11-87  23:50 

24 

25 

32 

6-14-87  19:33 _ 

19 

15 

16 

32 

33 

32 

KBRKEfS'RniflHIl 

19 

17 

22 

.  A  Va^sr 

Average  — - - 

^  *  TRVE 

0.96 

1.05 

RMS  Relative  Difl'erence 
- vs...  - 

0.09 

0.15 

estimate.  In  the  following  section  we  describe  an  approximate  solution  for  the 
dual  Gaussian  parameters  that  loads  to  an  efficient  velocity  estimation  algorithm. 


II.  Estimates  Using  /?(0)  and  R{r] 


The  number  of  unknowns  in  equations  (1)  can  be  reduced  from  eight  to  six  by 
assuming  that  the  ratio  of  high  and  low  beam  amplitudes  for  each  spectral  com¬ 
ponent  can  be  determined  from  a  known  beam  weighting  function, 


c  !  r\  “  1  f  l)"  1  , 

SiU)  =  wv,i— ==— exp  - — - 1  -t-  ic,-  o 

V27rcTi  2(7, 


(Z  0 


V^cr. 


-expl 


-U-hf 


2a‘ 


(2) 


Appendix  B  describes  a  solution  for  the  six  spectral  parameters  using  measure¬ 
ments  of  /?(0)  and  R{r)  in  the  two  receiving  beams  of  an  ASR.  The  important 
result  is  that  the  center  frequency  of  the  'low  altitude  induced"  spectral  com¬ 
ponent  can  be  derived  as: 


/?,(()) 

^2(0) 


Wi,2 

^*'2.2 


/?2(r)l 


(3) 


Thus  if  this  center  frequency  can  be  shown  to  accurately  represent  the  low  alti¬ 
tude  wind  field,  the  ASR  velocity  estimation  algorithm  could  be  simplified  to  a 
standard  pulse  pair  estimate,  preceded  by  linear  combination  of  the  low  and  high 
beam  autocorrelation  estimates.  This  rnetliod  was  in  fact  proposed  by  Weber  and 
Moser  (1987)  although  we  indicated  there  that  better  performance  would  be 
expected  if  the  high  and  low  beam  signals  wa're  "orthogonalized"  through  linear 
combination  of  their  in-phase  and  quadnUun'  components  prior  to  calculation  of 
autocorrelation  lags. 


As  described  in  Appendix  B,  the  [)recomputed  beam  weights  «',■  y  in  equation 
(3)  can  be  parameterized  by  an  upper  elevation  angle  (corresponding  to  the  height 
(jf  a  microburst  outilow)  ha-  the  low  altitude  spectral  component.  In  this  report  we 
will  treat  tliis  angle  as  range  independent,  although  belter  performance  miglit  be 
expected  were  it  a  decreasing  function  of  range,  h'or  the  inicrobursts  treated  ])re- 
viously,  'fable  II-2  comf)ares  velocity  shear  estimates  from  eciuation  (3)  with  those 
from  the  |)encil  bc-ain  weathei-  radar.  Upper  elevation  angles.  of  1".  2"  and  3" 
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were  used  in  computing  the  w^  j. 


Table  II-2:  Microburst 
_ ^ R(r] 

Velocity  Shear  Estimates  Using 

1  Dual  Gaussian  Fit _ _ _ 

Date/ Time 

9-10-87  22:47 

Pencil  Beam 
Radar  AV 
(m/s) 

23 

ASR  (^0  =  3°) 

13 

ASR  (00  =  2°) 

18 

ASR  (00  =  1°) 

26 

5-21-87  14:14 

28 

2B 

31 

36 

■  wm  I'l'i  iiT 

27 

22 

28 

■■1  m  rr  m'wm 

23 

18 

22 

27 

24 

17 

19 

25 

19 

12 

16 

22 

32 

33 

36 

39 

8-1-87  20:30 

19 

14 

18 

25 

Average  - - 

0.76 

0.91 

1.20 

RMS  Relative  Difference 

0.27 

0.15 

0.19 

Note  that  as  Oq  is  decreased,  the  weight  increases  (see  Appendix  B). 

When  wind  speed  decreases  with  altitude  as  is  normally  the  case  in  a  microburst, 
the  magnitude  of  the  phase  angle  of  the  high  beam’s  R{t)  autocorrelation  lag  will 
be  smaller  than  that  of  the  low  beam.  A  simple  vector  construction  illustrates 
that  these  observations  account  for  the  observed  inverse  relationship  between  Oq 
and  the  magnitude  of  the  velocity  shear  estimate  from  equation  (3). 

Since  the  heights  of  microburst  outflow  winds  vary  considerably  and  since  the 
above  events  were  at  different  ranges  from  the  radar,  the  "best"  value  for  was 
variable;  indeed  in  many  of  these  microbursts  the  best  velocity  measure  for  the 
approaching  core  corresponded  to  a  different  value  for  than  that  for  the  reced¬ 
ing  core.  Over  the  eight  events  however,  weights  y  computed  using  a  value  of 
2°  for  6q  resulted  in  a  shear  estimates  that  on  average  were  closest  to  those  meas¬ 
ured  by  the  pencil  beam  radar;  the  corresponding  RMS  relative  error  was  likewise 
minimum  for  this  setting.  While  the  match  to  the  pencil  beam  radar  measure¬ 
ments  could  have  been  improved  by  selecting  a  Oq  slightly  less  tlian  2.0,  we  felt 
that  such  fine  tuning  was  unwarranted  given  the  small  number  of  events  used  for 
the  evaluation. 

Comparison  of  Tables  II-l  and  II-2  suggests  that  the  simple  velocity  estimator 
of  equation  (3)  might  provide  comparable  accuracy  to  the  less  constrained 
approach  discussed  in  the  previous  subsection,  provided  that  a  suitable  weight 
wq  o/«'‘j,2  were  selected.  The  following  section  examines  the  extent  to  which  a  sin¬ 
gle  value  for  this  weight  would  provide  acceptable  velocity  estimates  for  a  much 
larger  data  set  than  evaluated  above. 
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ffl.  EVALUATION  OF  VELOCITY  ESTIMATES  USING  SIMULA¬ 
TIONS  AND  FIELD  MEASUREMENTS 


This  section  compares  ASR  radial  velocity  estimates  from  equation  (3)  to  those 
derived  from  the  full  spectral  differencing  technique  and  to  "truth"  as  defined  from 
microburst  models  (part  A)  or  simultaneous  measurements  from  the  pencil  beam 
weather  radar  (part  B).  The  data  processing  sequence  involves  the  elements 
described  in  Weber  (1987)  and  Weber  and  Noyes  (1988): 

(i)  clutter-map  based  high  pass  filtering  of  34  sample  sequences  from  the  high 
and  low  receiving  beams; 

(ii)  velocity  estimation  using  equation  (3)  or  the  spectral  differencing  algorithm 
described  by  Weber  and  Noyes  (1988).  For  the  autocorrelation  based 
method,  beam  weights  ir,  j  were  calculated  with  Oq  set  to  2.0"; 

(iii)  nine-point  nearest  neighbor  spatial  median  filtering  of  the  velocity  field  fol¬ 
lowed  by  smoothing  along  the  range  axis  with  a  five-point  Gaussian  filter; 

(iv)  divergence  detection  using  the  surface  outflow  portion  (Merritt,  1987)  of  the 
TDWR  microburst  detection  algorithm. 

All  data  used  in  part  B  were  collected  during  1988  in  wet  microbursts  near  Hunts¬ 
ville,  Alabama.  We  note  that  the  value  of  used  in  calculating  the  weighting 
coefficient  w'l.a/woo  was  derived  from  independent  microburst  data  collected  dur¬ 
ing  1987  as  described  in  the  preceding  section. 


A.  Velocity  Estimates  using  Simulated  ASR  Signals 

As  a  function  of  "true"  surface  outflow  radial  velocity.  Figures  111-1  compare 
the  bias  and  standard  deviation  of  velocity  estimates  from  the  spectral 
differencing  and  autocorrelation  based  algorithms  .  Estimate  standard  deviations 
are  shown  both  before  and  after  the  spatial  smoothing  described  above.  The  cal¬ 
culations  assume  that  the  outflow  velocity  is  constant  from  the  surface  to  100  m 
height;  the  radial  wind  then  changes  linearly  to  an  "upper  level"  velocity  that  is 
one-third  the  magnitude  of  the  surface  wind  and  opposite  in  direction. 
Reflectivity  factor  and  spectrum  width  are  taken  as  constant  in  altitude  witli 
values  of  10  dBz  and  2  m/s  respectively.  Figure  111-2  illustrates  this  mode!  for  a 
surface  outflow  velocity  of  15  m/s  towards  the  radar. 

The  velocity  estimate  performance  metrics  were  calculated  using  500  trials  of 
the  Monte  Carlo  signal  simulation  described  in  .\ppendix  C’.  Since  essentially  all 
data  from  our  experimental  ASR  in  Huntsville  were  obtained  at  a  uniform  pulse 
repetition  frequency,  the  simulations  here  assume  a  constant  PRF  of  980.s“'.  We 
show  in  Appendix  C  that  utilization  of  the  ASFC9’s  8/10  pulse  alternating  PRF 
waveform  would  produce  only  small  changes  in  velocity  estimate  accuracy.  Fig¬ 
ures  lll-l(a)  and  (b)  correspond  to  a  resolution  cell  where  the  weather  signal  to 
ground  clutter  ratio  is  sufficiently  large  that  high  pass  filtering  is  not  recjuired 
(Weber,  1987).  Part  (a)  assumes  that  the  rcsiolution  cell  of  interest  is  at  6  km 
range  and  [)art  (b)  considers  a  range  of  12  km. 

At  6  km  range,  biases  associated  with  either  method  are  less  than  1  m/s.  Esti¬ 
mate  standard  deviations  average  about  2  m/s  before  spatial  smoothing  and  are 
reduced  to  0.5  m/s  after  smoothing.  Note  that  our  spectral  differencing  algorithm 
(Weber  and  Noyes,  1988)  incoherently  averages  power  spectrum  density  (PSD) 
estimates  using  a  running  three  gate  average  in  range  before  calculating  radial 
velocity.  Thus,  more  extensive  spatial  smoothing  is  required  for  the  spectral 
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SPECTRAL  DIFFERENCING 


True  Velocity  (m/s) 


(a)  Range  of  6  km  and  no  clutter  filtering  are  assumed. 

Figure  IIl-l.  Spectral  differencing  and  autocorrelation  based  velocity  estimate  bias 
and  standard  deviation  versus  "true"  outflow  velocity.  The  velocity  model  of  F'ig- 
ure  III-2  is  scaled  proportionally  to  the  abscissae.  Solid,  dashed  and  chain  dashed 
curves  for  standard  deviation  pertain  respectively  to  single  resolution  cell  esti¬ 
mates,  estimates  after  spatial  median  filtering  and  after  Gaussian  smoothing  along 
the  range  axis. 


SPECTRAL  DIFFERENCING 


(c)  Range  of  6  kin  and  use  of  a  clutter  filter  are  a.ssvimed 


Figure  III-l.  (continued) 
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I'^igiirc  111-2.  Model  vortical  proliles  of  radial  velocity,  rollectivity  and  spectrum 
width  in  a[)proaching  |)ortion  of  a  microburst. 


differencing  algorithm  than  for  the  autocorrelation  based  method  to  achieve  the 
stated  estimate  variance. 

Near  zero  Doppler,  the  standard  deviation  of  the  spectral  differencing  estimates 
increases  slightly.  In  our  model,  there  is  no  v'^ertical  shear  in  the  radial  wind  for 
zero  velocity  (i.e.  at  the  microburst  center).  Without  this  shear,  underlying  high 
and  low  beam  PSDs  are  the  same  and  the  placement  in  the  velocitv  domain  of  the 
low-high  beam  difference  spectrum’s  positive  lobe  is  governed  by  statistical 
fluctuations  in  the  spectral  estimates.  This  lobe  —  which  is  used  to  generate  the 
low  altitude  velocity  estimate  —  can  occur  anywhere  where  the  PSDs  are  within  a 
few  decibels  of  their  maximum  values.  This  circumstance  accounts  for  the 
increase  in  estimate  standard  deviation  at  low  Doppler.  Cionversely,  for  the  auto¬ 
correlation  based  algorithm,  velocity  estimate  standard  deviations  are  minimum 
at  zero  Doppler  where  the  integrated  spectrum  width  is  smallest.  This  is  con¬ 
sistent  with  Zrnic's  (1977)  findings  that  the  conventional  weather  radar  pulse-pair 
mean  velocity  estimate  variance  decreases  with  decreasing  spectrum  width. 

The  larger  biases  at  12  km  range  (Ill-l  (b))  represent  an  underestimate  of  the 
surface  w'ind  s{)ee(i:  the  error  increases  as  surface  velocity  and  vertical  shear 
become  larger  in  magnitude.  Maximum  calculated  bias  for  the  spectral 
dift'erencing  and  autocorrelation  l>ased  algoritl)ms  are  3  m/s  and  1  m/s  respec¬ 
tively.  Standard  deviations  are  slightly  larger  than  at  5  km  range. 

In  Figure  111- 1(c)  we  repeat  the  6  km  range  calculation  assuming  that  high  pass 
clutter  filtering  as  described  by  Anderson  (1987)  or  Weber  (1987)  is  reciuired.  The 
filter's  stop  band  width  is  ±-1.7  m/s.  providijig  39  dB  allenuation  of  scan  niodn- 
lated  ground  clutter.  In  addition  to  removing  l<')vv  Doppler  signal  |)ower.  the  17- 
coefl'icient  filter  reduces  the  number  of  valid  data  points  for  velocity  estimation 
from  3f  to  18.  We  did  not  include  a  ground  clutter-induced  spectrum  component 
in  our  signal  simulation,  'bhus.  simulated  weather  spectrum  distortion  may  be 
larger  than  would  occur  with  real  data  wlu're  the  additiot)al  low  Do|)pler  power 
from  grotind  clutter  would  at  least  partially  offset  the  filtering. 


'Fhe  overestimate  of  velocity  magnitude  resulting  from  spectrum  distortion  is 
maximum  when  the  surface  velocity  approaches  the  e(lges  of  the  lilter  stoj)  band: 
at  this  velocity  bia.ses  are  about  2  m/s  for  either  algorithm.  Fstimate  standard 
ileviations  are  also  large  at,  low  Dojipler  where  most  of  the  weather  i)owcr  sp(  c- 
trum  lies  inside  the  filter's  slop  b:ind;  evi-n  after  sjiatial  smoothing,  estimate  stan¬ 
dard  deviation  would  be  about  2  m/s  for  low  Dopjiler  signals.  At  higher  Dopider 
velocities.  the  effect  of  the  smalh'r  number  of  available  data  sam|)les  is  evident  in 


increa-sed  estimate  variance  relative  to  th(‘  no-lilter  sit uat ion; 


incix'ase  is  mori' 


pronounceil  for  the  autocorrelation  based  algorithm.  Outside  the  low  Doppler 
interval  however,  spatial  smoothing  is  effective  in  redm’ing  the  standard  deviation 
for  either  method  to  less  than  I  m/s. 


Figure  1II-.3  plots  velocity  estimate  bias  and  standard  deviation  as  a  function  of 
microbui’st  reflectivity  factor,  d'he  calculation  us(‘d  the  profile  of  liguia'  111-2  (sur¬ 
face  velocity  -13  m/s)  excejit  that  the  n-lleetivity  factor  was  varied  from  30  dBz 
down  to  -10  dBz.  We  assumed  a  |•ange  of  (i  km  and  a  sensitivity  time  control 
function  that  jilaces  the  receiver  noise  level  at  an  <‘<|uivalent  weather  reflectivity 
factor  of  0  dBz  (.s('e  Weber  and  Moser.  1987).  When  tlu'  reflectivity  factor  is 
greater  than  0  dBz.  l)ot  h  the  spectral  dilfer<>ncing  and  autocorrelation  liased  velo¬ 
city  ('stimates  exhibit  minimal  bias  and  standard  deviations  less  than  1  m/s  after 
siiatial  smoothing.  .\s  the  signal  to  noise  ratio  Ix-comes  negative,  bias  and 
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standard  deviation  increase  rapidly.  The  positive  bias  reflects  the  trend  of  the 
velocity  estimates  towards  0  m/s  when  noise  dominates  the  spectrum.  Note  that 
this  calculation  does  not  consider  the  impact  of  ground  clutter.  Weber  (1987) 
showed  that  ground  clutter  may  prevent  accurate  velocity  measurements  at  short 
range  when  the  reflectivity  factor  is  less  than  20  dBz. 

The  dependence  of  velocity  estimate  accuracy  on  the  vertical  extent  of  the 
outflow  is  illustrated  in  Figure  III-4.  Here,  bias  and  standard  deviation  are  plot¬ 
ted  as  functions  of  outflow  height,  assuming  a  range  of  6  km.  The  radial  velocity 
profile  is  that  of  III-2  scaled  in  height  by  a  factor  varying  from  0.25  to  2.25.  The 
shallowest  outflow  considered,  therefore,  has  maximum  radial  winds  extending  25 
m  above  the  surface  and  reaches  the  "up|)er  level"  velocity  at  250  ni.  For  this 
height,  both  algorithms  result  in  significant  underestimates  of  the  near  surface 
radial  wind  magnitude.  The  bias  is  fi  and  7  m/s  respectively  for  the  spectral 
differencing  and  autocorrelation  based  algorithms.  With  increasing  outflow  depth, 
the  bias  decreases  rapidly  to  valiums  less  than  1  m/s  when  the  outflow  layer  is  75 
m  deep.  Estimate  standard  deviation  increases  by  about  .50  percent  for  very  shal¬ 
low  outflows. 

I'or  Huntsville  microbursts  in  1987,  Weber  and  Noyes  (1988)  plotted  the  distri¬ 
bution  of  heights  at  which  the  radial  v^elocity  dropped  to  half  of  its  maximum 
value.  The  median  value  for  this  "half-height"  was  3.50  m  and  about  15  percent  of 
the  microbursts  exhibited  half-heights  less  than  200  m.  These  data  were  com¬ 
piled  from  events  centered  as  far  as  12  km  from  the  radar  where  the  3  dB  beam 
width  spanned  300  m.  Thus,  angular  resolution  may  have  been  inadequate  for 
measuring  the  actual  depths  of  some  of  the  microbursts.  Biron  and  Isaminger 
(1989)  analyzed  RHI  scans  of  microbursts  wit  hin  8  km  of  Lincoln  Ijaboratory’s  H 
pencil  beam  weather  radar.  Vertical  resolution  was  therefore  140  m  or  better.  The 
median  half-height  for  microbursts  they  measured  during  1986  in  Huntsville  was 
400  m;  all  events  were  between  300  and  1000  m  deep  by  this  measure.  For 
Denver  microbursts  observed  during  1987,  the  half-height  distribution  extended 
from  200  to  1100  m  with  a  median  value  of  600  m.  The  median  height  of  max¬ 
imum  velocity  in  both  locales  was  within  the  lowest  200  m  .ACL.  For  comparison 
with  these  statistics,  the  upper  abscissa  labels  on  111-4  give  the  half-height  for  the 
vertical  wind  profile  wc  assumed  in  our  calculations.  At  400  m  half-height,  esti¬ 
mate  bias  is  less  than  I  m/s  for  either  method:  the  shallowest  15  ])ercent  of 
(Jilt flows  measured  in  Huntsville  during  1987  would  be  subject  to  biasc's  2.5  m/s  or 
greater  at  6  km  range. 

Radial  velocity  versus  range  signatures  through  a  model  microburst  are  simu- 
lateil  in  i''igur('s  111-5.  The  basic  prolih'  is  again  that  of  111-2  but  the  surface 
(Jilt  flow  velocity  vvas  varied  siiiusiddally  with  range  and  upper  level  winds  were 
scab'll  pr(jp(jrtionally.  Thus  the  (jverall  structure  exhibits  surface  divergence  with 
comp('nsating  c(jnvergence  abjft  as  is  characi eristic  of  measured  microburst  wind 
fields.  The  moders  maximum  approaching  and  reci'din^  vebjcity  coits  are 
se[)araled  by  2  km.  In  each  plot.  tlies(jfid  line  is  (In'  "l  rue  '  surface  radial  wind 
pattern:  d.ashed  and  chain-dashed  curve's  are' single  re'alizations  of  the  velejcity  sig¬ 
nature  estimated  with  the  s()e’ctral  eliflere>iie-ing  and  autejeorrelat ion  based  algo¬ 
rithms. 

In  lll-.5(a)  and  (b)  the  "microburst"  is  e-e'iitered  at  6  and  12  km  respect ivi'ly:  the' 
iijiper  panels  in  each  figure'  simulate  vebjcity  e'stimatt's  wlu're  no  spatial  snujothing 
has  been  ajipfied.  'I'lie  variance  of  the'  iinsnujot  he'd  ve’lejcity  e-stiniati's  eiisru|)ts 
t  h('  moiKjtonically  increasing  patte'rn  bi'lwe'en  the'  velocity  extrema  to  the  e'Xti'ut 
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Figure  I1I-4.  Spectral  differencing  and  autocorrelation  bailed  velocity  estimate 
bias  and  standard  deviation  versus  outflow  height.  The  velocity  model  of  Figure 
III-2  is  scaled  proportionally  to  the  lower  abscissa  labels.  The  curve  for  standard 
deviation  pertains  to  estimates  after  spatial  median  filtering  followed  by  Ciaiissian 
smoothing  along  the  range  axis. 
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(a)  Microburat  is  assumed  lo  be  centered  at  6  km  range. 

Figure  1II-5.  Simulated  microl)urst  radial  velocity  versus  range  mea.suroinents 
from  an  ASR.  The  solid  sinusoid  is  the  model  for  the  surface  radial  wind  associ¬ 
ated  with  the  microburst.  Dashed  and  chain  dashed  curves  represent  estimates 
using  the  spectral  differencing  and  autocorrelation  based  algoritbms.  'I'lie  upjicr 
plot  simulates  single  resolution  cell  estimates;  the  lov/er  plot  ineludes  the  spatial 
filtering  described  in  the  tt'xt. 
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(b)  Microhiirst  is  assumed  (o  be  eeiilered  at  12  km  range 
Figure  III-5.  (continued) 
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that  the  TDWR  microburst  detection  algorithm  might  not  identify  shear  seg¬ 
ments.  The  lower  panels  simulate  the  spatial  smoothing  described  above.  The 
smoothed  velocity  fields  correlate  well  with  the  model  wind  pattern  although  the 
bias  indicated  in  Figure  Ill-l(b)  is  evident  when  the  microburst  is  centered  at  12 
km  range. 

This  simulation  based  analysis  indicates  that  velocity  estimates  from  the  auto¬ 
correlation  based  and  spectral  differencing  algorithms  exhibit  very  similar  bias  and 
standard  deviation.  A  representative  value  for  standard  deviation  after  spatial 
filtering  is  0.5  m/s.  The  most  significant  factor  determining  the  bias  associated 
with  either  algorithm  is  the  rate  of  change  of  radial  velocity  with  elevation  angle. 
For  very  shallow  microbursts  ("half-height"  less  than  200  m)  and  microbursts 
beyond  about  12  km  range,  differential  velocity  underestimates  of  15  percent  or 
more  could  be  expected.  However,  for  the  majority  of  microbursts  in  the  opera¬ 
tionally  significant  region  within  12  km  of  an  ASR,  our  analysis  indicates  that  the 
accuracy  of  either  velocity  estimator  is  good. 
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B.  Field  Measurements  from  Huntsville  Experimental  ASR 

Lincoln  Laboratory’s  airport  surveillance  radar  weather  detection  experiment  in 
Huntsville,  Alabama  during  1987  was  described  by  Weber  and  Noyes  (1988). 
Equipment  and  operating  procedures  during  the  summer  of  1988  were  identical 
except  that  the  scan  rate  of  the  C-band  pencil  beam  radar  used  for  "truth"  was 
increased.  The  basic  scan  pattern  for  microbursts  in  1988  consisted  of  two  360® 
PPI  scans  at  elevation  angles  of  0.6®  and  1.5®,  followed  by  two  RHI  scans  through 
the  outflow.  This  sequence  was  repeated  at  one  minute  intervals. 

Figures  III-6  compares  images  of  the  radial  velocity  field  estimated  from  our 
airport  surveillance  radar’s  signals  with  that  measured  by  the  pencil  beam 
weather  radar.  Data  arc  from  a  microburet  producing  thunderstorm  on  15 
August  1988.  The  upper  left  panel  is  the  pencil  beam  measurement  from  a  scan 
at  0.6®  elevation  angle.  Two  microbursts  were  present,  a  strong  outflow  centered 
at  10  km  range/130®  azimuth  and  a  weaker  event  at  15  km/65®.  ASR  estimates 
using  the  low-high  beam  spectral  differencing  technique  are  shown  in  the  upper 
right  panel  with  the  corresponding  autocorrelation  based  estimate  in  the  lower 
left.  High  and  low  beam  signals  from  the  same  antenna  scan  were  employed  for 
these  estimates.  The  autocorrelation  based  estimate  obtained  when  the  high  beam 
signal  was  collected  from  the  following  antenna  scan  is  shown  in  the  lower  right. 
For  current  ASRs,  this  mode  of  data  collection  would  be  required  for  operation 
using  circular  polarization. 

The  velocity  fields  derived  from  the  ASR  signals  are  in  good  agreement  with 
that  measured  by  the  weather  radar.  In  particular  the  presence  of  the  two  micro¬ 
bursts  is  clearly  indicated  and  the  ASR  velocity  differential  estimates  are  within  1 
m/s  of  the  pencil  beam  measurements.  The  spectral  domain  and  autocorrelation 
based  wind  field  estimates  from  the  ASR  signals  are  likewise  in  good  agreement, 
particularly  within  the  microbursts.  Collection  of  high  and  low  beam  autocorrela¬ 
tion  estimates  on  alternate  antenna  scans  (lower  right  panel)  did  not  significantly 
change  the  velocity  estimates. 

We  suspect  that  the  1.4®  beamwidth  of  the  weather  radar  was  too  large  for 
accurate  measurement  of  the  receding  outflow  component  in  the  microburst  to  the 
southeast.  As  in  .some  of  the  other  cases  presented  below,  the  stronger,  more 
homogeneous  receding  volume  depicted  in  the  ASR-based  fields  may  well  be  a 
better  representation  of  the  actual  wind  field. 

Additional  examples  are  shown  in  F'igures  II1-7  through  111-9  .  In  each  figure, 
the  upper  panel  displays  the  pencil  beam  radar  0.6®  velocity  field.  The  lower  left 
and  lower  right  panels  are  ASR-based  estimates  using  respectively  the  autocorrela¬ 
tion  and  spectral  differencing  approaches.  High  and  low  beam  signals  were  from 
the  same  antenna  scan  since  we  do  not  normally  transfer  data  from  adjacent  scans 
for  analysis. 

These  examples  again  indicate  good  agreement  amongst  the  pencil  beam  and 
ASR  derived  velocity  fields.  For  the  displayed  scans  on  21  June  and  25  June, 
both  ASR  based  velocity  difl'erential  estimates  are  within  2  m/s  of  the  pencil  beam 
radar’s  measurement.  The  velocity  ‘^liear  estimate  from  the  autocorrelation  based 
ASR  field  de[)icted  in  Figure  III-8  is  6  m/s  (23^o)  larger  than  that  derived  from 
the  |)encil  beam  radar;  the  spectral  differencing  approach  resulted  in  a  3  m/s 
"overestimate"  in  this  case.  Note  that  the  the  divergent  outflows  shown  in  these 
examples  extend  to  ranges  as  large  as  20  km  (I'figures  111-8  and  III-9). 
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I'igiirc  HI-7.  Images  of  the  radial  velocity  field  in  a  microburst  producing  thunderstorm  near  Huntsville,  Ala- 
I  :iiu:i  on  21  June  1988.  The  upper  panel  show.s  measurements  from  the  MIT  C-band  beam  weather  radar 
caniiiiig  ill  PIT  mode  at  OG"  elevation  angle.  Ixjwer  right  panel  is  field  estimated  from  ASR  signals  using 
I  he  spfctral  dilTenuicing  method.  Lower  left  panel  is  corresponding  field  from  autocorrelation  based  algo 
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igtire  II1-8.  Images  of  the  radial  velocity  field  in  a  microburst  producing  thunderstorm  near  Huntsville,  Ala- 
.'ima  on  24  June  1988.  The  upper  panel  .shows  measurements  from  the  MIT  C-band  beam  weather  radar 
aiuiing  in  PPI  mode  at  0.6°  elevation  angle.  Lower  right  panel  is  field  estimated  from  ASR  .signals  using 
Ih'  spectral  dilTerencing  method.  lx)wer  left  j)anel  is  corresponding  field  from  autocorrelation  based  algo 
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I’igure  111-9.  Images  of  radial  velocity  field  in  a  microburst  producing  thunderstorm  near  Huntsville,  Alabama 
III  ‘25  June  1988.  The  upper  panel  shows  measurements  from  the  MIT  C-band  beam  weather  radar  scanning 
in  Pl’I  mode  at  0  6°  elevation  angle.  Lower  right  panel  is  field  estimated  from  ASR  signals  using  the  spectral 
lilh  ri  iM'ing  method.  Lower  left  panel  is  corresj)onding  field  from  autocorrelation  based  algorithm. 


A  more  extensive  evaluation  of  the  accuracy  of  the  ASR  velocity  estimates  was 
derived  by  "scoring"  detections  from  the  TDWR  surface  outflow  detection  algo¬ 
rithm  against  truth  as  determined  from  manual  examination  of  the  pencil  beam 
weather  radar  data.  The  procedure  and  associated  scoring  rules  are  described  in 
Weber  and  Noyes  (1988).  Our  evaluation  utilized  approximately  600  scans  from 
the  airport  surveillance  radar  in  1988,  taken  during  35  microbursts  on  13  separate 
days.  The  scoring  was  confined  to  microbursts  centered  within  the  operationally 
significant  area  extending  12  km  from  the  radar. 

Table  III-l  summarizes  the  results  of  scoring  on  a  scan  by  scan  basis.  The 
listed  performance  metrics  are: 

(i)  probability  of  detection  —  the  number  of  detected  rnicroburst  signatures 
divided  by  the  total  number  of  microburst  signatures: 

(ii)  probability  of  false  alarm  —  the  number  of  algorithm  alarms  not  associated 
with  microbursts  divided  by  the  total  number  of  alarms; 

(iii)  bias  ~  the  average  difference  between  ASR-based  and  pencil  beam  radar 
microburst  differential  velocity  estimates.  This  is  expressed  both  in  absolute 
units  (ni/s)  and  relative  to  the  pencil  beam  radar  Al'/j  measurement; 

(iv)  root  mean  squared  (RMS)  difference  between  the  pencil  beam  radar  and 
ASR-based  velocity  differential  estimates.  This  equals  the  square  root  of  the 
sum  of  squared  estimate  bias  and  estimate  variance. 

These  metrics  are  tabulated  separately  for  all  microbursts  and  for  microbursts 
with  differential  velocities  greater  than  15  and  20  m/s. 


The  results  confirm  the  favorable  prognosis  for  an  ASH's  capability  to  detect 
wet  microbursts  that  we  derived  from  analysis  of  data  collected  during  1987 
(Weber  and  Noyes.  1988).  Both  velocity  estimation  algorithms  supported  detec¬ 
tion  and  false  alarm  probabilities  within  the  0.9/0. 1  bounds  called  for  by  the  I'AA 
in  its  TDWR  system  recpiirements  statement.  Over  all  microbursts  seoi’cd.  tlu' 
dual-beam  autocorrelation  metho(J  I’csiilted  in  somewhat  reduced  (h'tection  jx'rfoi'- 
mance  relative  to  full  spectral  differencing:  (h-tc'clion  probability  decreased  by  two 
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percent  while  the  false  alarm  probability  and  discrepancy  with  pencil  beam  radar 
AVff  measurements  increased.  These  differences  disappeared  when  scoring  wa.s 
restricted  to  the  subset  of  microburst  scans  where  the  velocity  differential  exceeded 
20  m/s. 

Detailed  analysis  of  these  data  indicates  that,  for  either  velocity  estimation 
method,  eighty  percent  of  the  missed  detections  were  associated  with  recognizable 
(by  a  human  observer)  divergence  patterns  in  the  ASR-based  radial  velocity  field. 
Tbe  automatic  microburst  detection  algorithm  did  not  declare  an  alarm  in  these 
cases  because  its  spatial  size  and/or  continuity  requirements  were  not  met  (Mer¬ 
ritt,  1987).  The  remaining  misses  occurred  in  weak  microbursts  (AT/;  <  15m /.s ) 
where  the  ASR  based  velocity  differential  estimates  did  not  exceed  the  10  m/s 
microburst  threshold.  Our  analysis  of  missed  detections  does  not  suggest  that 
divergent  outflows  are  depicted  more  clearly  in  spectral  differencing  versus  auto¬ 
correlation  based  velocity  fields.  Given  the  subtleties  of  the  automatic  microburst 
detection  algorithm,  we  regard  the  2  percent  overall  difference  in  detection  proba¬ 
bilities  using  the  two  velocity  estimators  as  insignificant. 

Almost  all  false  alarms  occurred  in  regions  where  the  pencil  beam  radar  meas¬ 
ured  divergence,  but  below  the  required  threshold.  As  seen  from  Table  III-l,  auto 
correlation  based  estimates  of  velocity  divergence  were  larger  on  average  than 
measured  by  the  pencil  beam  radar,  particularly  for  weak  (AV'yj  15m /.s )  micro¬ 
bursts.  This  "bias"  produced  the  higher  overall  false  alarm  rate  associated  with 
that  estimator. 

Our  simulations  (Section  III-A)  predict  that  biases,  if  present,  should 
correspond  to  underestimates  of  velocity  differential  and  should  be  larger  for 
strong  microbursts  where  the  vertical  gradient  in  radial  wind  speed  is  highest. 
This  is  clearly  inconsistent  with  Table  lll-l.  A  possible  exi)lanation,  as  alluded  to 
previously,  is  that  the  beamwidth  and/or  ground  clutter  suppression  capability  of 
the  pencil  beam  radar  were  inadequate  for  accurate  measurement  of  (he  strongest, 
near  surface  microburst  winds.  A  uniform  increase  in  "true"  velocity  differential 
for  the  1988  data  set  would  at  least  produce  the  exi)ected  trend  for  bias  versus 
inicroburst  intensity. f 

Overall,  these  statistics  confirin  the  previous  simulations  and  case  studies  indi¬ 
cating  that  the  dual-beam  autocorrelation  based  velocity  estimate  of  equation  (3) 
should  support  wet  microburst  detection  at  a])|)roximately  the  same  level  of 
confidence  as  would  be  obtained  through  full  spectral  differencing.  Examination 
of  the  storm  cases  used  for  the  detection  performance  statistics  in  Table  lll-l  is 
continuing.  The  analysis  will  f)rovide  more  detailed  understanding  of  tnissed 
detections  or  false  alarms,  and  the  small  dilferences  in  detection  algorithm  |)erfor- 
niance  observed  using  the  two  velocity  estimators. 


^  Mca.siire(l  bia.s  for  the  s|)(Ttral  (lifTcrfririnK  algoritlim  ai)|)li('cl  lo  our  1987  data  sot  was 
qiialitativcly  ronsi.stcnt  witli  simulations  in  Section  III-A.  Compare  for  e.xami'le  Figures 
lll-l(h)  from  this  report  and  the  "shear  ratio"  (i.e.  l)ias)  plotted  in  Figure  \’I-15  of  Weher 
and  Noyes  (1988), 
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IV.  DISCUSSION  AND  FUTURE  WORK 


The  analysis  herein  confirms  our  previous  assessment  that  a  suitably  modified 
airport  surveillance  radar  would  provide  an  operationally  useful  capability  for 
automatic  detection  of  "wet"  microbursts.  In  this  report,  we  developed  an 
efficient,  autocorrelation  based  low  altitude  velocity  estimator  based  on  the 
assumption  that  the  pow'er  spectrum  of  weather  echoes  measured  by  an  ASR 
could  be  adequately  represented  with  a  bimodal  (Jaussian  model.  Evaluation  of 
microburst  detection  algorithm  performance  using  tlie  resulting  velocity  fields 
indicated  that  accuracy  comparable  to  the  most  effective  velocity  estimator 
treated  by  Weber  and  Noyes  (1988)  —  low-high  beam  spectral  differencing  —  might 
be  achieved  with  considerable  reduction  in  computational  and  hardware  require¬ 
ments.  The  combination  of  high  confidence  and  computational  efficiency  make 
this  technique  an  attractive  initial  candidate  for  implementation  in  our  real  time 
wind  shear  processor. 

Ongoing  evaluation  of  ASR  velocity  ('sUmation  techniques  will  seek  to  further 
quantify  tradeoffs  between  accuracy  and  processing  complexity.  In  addition  to  the 
autocorrelation  based  and  spectral  differencing  methods  treated  here,  we  are  con¬ 
tinuing  to  assess: 

(i)  low  beam  only  estimators,  for  example  pulse  pair  processing  following  high 
pass  filtering.  Necessary  front  end  modifications  would  be  reduced  if  an  .\SR 
wind  shear  processor  only  accessed  signals  from  the  low  receiving  beam: 

(ii)  mapping  from  the  phase  of  the  cross  spectral  density  between  low  and  high 
beam  signals  to  elevation  angle  (Anderson,  1989).  This  teehnique  may  pro¬ 
duce  more  accurate  near-surface  velocity  estimates  and  offers  the  potential  for 
three-dimensional  rellectivity  and  wind  measurements.  Note  however  that 
high  and  low  beam  signals  must  be  accessed  simultaneously:  thus  on  current 
.■\Sfls  this  method  could  not  be  employed  during  o])eration  with  circular 
polarization. 

Each  of  the  candidate  algorithms  will  be  scored  against  all  available  tiata  from 
our  field  experiments  in  Huntsville  and  our  current  site  near  Kansas  City,  Mis¬ 
souri.  Our  goal  is  to  develop  reliable  performance  statistics  in  order  to  specify  the 
design  of  an  ASR  wind  shear  processor. 

We  did  not  consider  here  the  use  of  rellectivity  factor  measurements  from  an 
ASR  as  a  sui)porting  and/or  precursory  indicator  of  microburst  activity.  Eield 
measurements  with  pencil  beam  I)o[)[)ler  weather  radars  have  shown  that  descend¬ 
ing  reflectivity  cores  fre(|uently  pixsage  the  develo[)mcnt  of  strong  stirface  outflows 
(Isaminger,  1988).  While  temporal  growth  i)f  the  rellectivity  field  measured  by  the 
low  beam  of  an  .'\S|-1  may  fu-ovide  indirect  evidence  of  a  descending  core,  less 
ambiguous  information  could  be  be  derived  from  cotnparison  of  rellectivity  in 
"upper"  and  "lower"  beams,  d’hf'  high  and  h)w  beams  could  be  used  directly  in 
this  matiner  although  they  overlap  substantially.  Better  differentiation  would  lie 
provided  by  combining  clata  from  the  two  beams  in  an  effort  to  explicitly  separate 
received  power  itito  upper  and  lower  altitmh'  components,  d'his  could  be  done  in 
the  freciuency  domain  using  techniciues  analogous  to  those  described  for  velocity 
estimation.  Alternately,  amplitudes  from  the  dual  (iaussian  si)ectral  model  —  cal¬ 
culated  as  in  the  appetidices  to  this  report  —  could  provide  the  desired  relh'ctivity 
measures.  We  an-  examining  the  ihwelopment  of  the  reflectivity  field  as  seen  by 
our  exp(>riment al  ASR  to  determine  ai)pro|)riate  data  processing  a|)|)roaches  for  its 
use  in  microburst  detection. 


As  pointed  out  by  Weber  and  Moser  (1987),  slightly  lower  gain  and  increased 
"beamfilling  loss"  for  an  ASR’s  high  beam  reduce  sensitivity  to  low  reflectivity 
weather.  This  may  limit  the  utilization  of  dual  beam  techniques  for  "dry"  micro¬ 
burst  and  gust  front  detection.  As  currently  implemented,  our  dual-beam  velocity 
estimation  algorithms  compare  received  power  from  both  beams  to  receiver  noise; 
if  only  the  low  beam  signal  passes  this  threshold  test,  we  revert  to  a  standard  lo\.- 
beam  mean  velocity  estimate  for  that  resolution  cell.  Analyses  of  gust  fronts 
observed  in  Huntsville  and  simulation  of  dry  microbursts  using  volume-scan  pencil 
beam  weather  radar  data  are  underway  to  quantify  the  ability  of  ASRs  to  detect 
low  reflectivity  wind  shear.  These  will  establish  firm  low-end  reflectivity  limits  for 
ASR  wind  measurements. 

We  noted  previously  instances  where  the  accuracy  of  our  pencil  beam  "truth" 
radar’s  velocity  field  was  questioned  owing  to  its  relatively  broad  1.4"  beam.  Lim¬ 
ited  clutter  suppression  capability  and  slow  scanning  also  reduce  confidence  in 
velocity  measurements  from  this  radar.  To  more  reliably  quantify  the  accuracy  of 
ASR  wind  estimation  algorithms,  we  have  located  our  experimental  ASR  at  the 
same  site  as  Lincoln  Laboratory’s  TDWR  test  bed.  That  system  will  provide  a 
narrower  beam  (l®  in  1989  and  0.5®  after  conversion  to  C-band  operation  in 
1990),  good  clutter  suppression  and  rapid  volumetric  scanning.  The  scan  strategy 
facilitates  direct  comparisons  with  an  ASR’s  surface  wind  estimates  and  under¬ 
standing  of  errors  through  analysis  of  the  three-dimensional  wind  field.  Field 
measurements  will  continue  in  1990  at  Orlando,  Florida.  Collection  and  analysis 
of  ASR  wind  measurements  in  these  varied  environments  will  refine  our  under¬ 
standing  of  the  capabilities  of  ASRs  for  wind  shear  detection. 
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APPENDIX  A:  ESTIMATION  OF  DUAL  GAUSSIAN  SPECTRAL 
PARAMETERS  USING  R{t)  and  R{2t) 

The  first  three  autocorrelation  lags  associated  with  the  assumed  dual  Gaussian 
spectra  in  equation  (1)  are: 

/?,(0)  =  a.-.i+fli  .2  (^'1) 

/?,(r)  =  a,-  iexp(— 27rVfT^)exp(2  27r/  ,r)  +  a,- 2exp(— 27r“(72“7’“)exp(2  27r/o7')  (A-2) 

/?,{2t)  —  a,-  iexp(— 87rVfr)exp(*47r/  ,r)  +  a,_2exp(— 87i^o'2r^)exp(2  47r/2r)  (A-3) 

where  i=l,2  indexes  the  low  and  high  beams. 


Measurement  of  /?,(0)  and  /fi(r)  provide  only  six  relations  for  the  eight  unk¬ 
nown  spectral  parameters  whereas  the  addition  of  /?,(27)  results  in  an  over¬ 
determined  system  of  equations.  We  derive  therefore  a  solution  based  on  measure¬ 
ments  of  the  autocorrelation  function  at  lags  t  and  2r.  We  first  eliminate 
between  the  two  beams  in  equations  (A-2)  and  (A-3)  the  autocorrelation  function 
components  due  to  one  of  the  spectral  modes.  For  example: 


R,[r)-^RiT) 
a-i  o 


jexp(— 2;rV,^r-)exp(z'27r/  ,7-) 


(A-4) 


/?i(2r) - ^/?2(2r)  = 

d  O  O 


^  l,lf^2,2~“  1.2'* 2,1 


'exp(— 87rVfr“)exp(i  47r/  jf)  (A-5) 


The  center  frequency  /,  of  this  spectral  component  may  be  determined  itera¬ 
tively  by  finding  the  weight  (ai,2/«2.2)  ^hat  the  phase  angle  of  equation  (A-5) 
is  twice  that  of  (A-4).  The  phase  angle  of  equation  (A-5)  may  have  to  be  unfolded 
before  comparison  to  that  in  (A-4).  In  performing  the  .search,  we  make  use  of  the 
known  high  and  low  beam  patterns  to  constrain  the  possible  values  of  this  weight. 
Our  convention  is  that  the  first  Gaussian  component  in  equation  (1)  is  associated 
with  low  altitude  winds  (below  an  elevation  angle  9q)  and  the  second  with  winds 
aloft.  Thus  the  desired  weight: 


fZ{0)Bt{0)Bi{e)d9 

On 


JZ{9)B,{9)B2{0)d9 

^0 


can  be  shown  to  be  in  the  range  0  to  1.0  for  reasonable  values  of  9^)  and  the 
reflectivity  factor  profile  Z{9).  Here  B^  and  By  are  the  low  and  high  beam  one¬ 
way  elevation  antenna  patterns. 


Having  determined  the  weight  ((ii,2/«2.2)  'Iip  spectrum  width  a,  is  calculated  by 
dividing  the  magnitude  of  (A-4)  by  that  of  (A-5): 


|/?,(r)-^/^,(r)| 


(A-fi) 
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The  factor  <^i.2Q2,i|  readily  derived  from  either  (A-4)  or  (A-5). 

«2,2 

An  analogous  procedure  is  used  to  derive  the  center  frequency  and  width  of  the 
"upper  level"  Gaussian  component  as  well  as  the  ratios  (01,1/02,1) 

[  ”i.i*^2,2-«i^,2«2_£|^  functions  of  the  a,j  have  now  been  determined;  from 

®2,1 

these  the  values  of  these  am.plitudes  can  be  readily  determined. 


42 


APPENDIX  B:  ESTIMATION  OF  DUAL  GAUSSIAN  SPECTRAL 
PARAMETERS  USING  R{0)  and  R{t) 


We  now  assume  that  the  two  spectral  components  have  intrinsic  amplitudes  ay 
that  are  modified  by  known  beam  weighting  coefficients  w.  y.  The  weighting 
coefficients  are  determined  by  equating  the  integrated  power  of  the  spectrum  in 
equation  (2)  with  that  which  would  be  determined  from  a  known  profile  of  the 
radar  reflectivity  factor  Z{6): 

Ri{0)  =  CiW,-  1+0211',- 0 

TT 

y 

=  jz{0)B^{e)Bi(e)dd  (B-i) 

0 

n 

^0  T 

=  z,_,,jB,{e)B,{9)de  +  Zo_,,jB,{e)B,{9)de 
0  00 


Note  that  the  radar  constant  is  taken  as  unity  in  equation  (B-1).  The  amplitudes 
ay  are  seen  to  be  "average"  reflectivities  within  assumed  lower  (below  ^o)  ^rid 
upper  elevation  angle  intervals.  The  y  are  integrals  between  the  appropriate 
elevation  angle  limits  of  the  low  or  high  beam  two-way  elevation  patterns. 

The  first  two  autocorrelation  lags  of  the  the  spectra  in  equation  (2)  are  given  by 
the  first  line  in  (B-1)  and: 

/?,(r)  =  iaiexp(— 27r“cri-r)exp(2  27r/ ,r)  +  2fl2®^P(~27rV2r)exp(?  27r/.2r)  (B-2) 

The  six  unknown  parameteis  are  now  the  ay,  cry  and  /y  of  the  two  Gaussian  spec¬ 
tral  components. 


The  amplitudes  ay  are  easily  computed  by  applying  the  Inverse  of  the  precom¬ 
puted  beam  weighting  matrix  tc,  y  to  the  measured  /f,  (0).  Center  frequencies  and 
widths  of  the  two  spectral  components  can  be  determined  by  eliminating  the  other 
component  between  the  high  and  low  beams  in  equation  (B-2).  For  example: 

R,ir)-:^R,{r)  =  |  ,^,^xp(-2;rVfr)exp(»2;r/ir)  (B-3) 


Wo  •) 


Woo 


The  pha.se  angle  of  (B-3)  determines  the  center  frequency  /  ,.  The  spectrum  width 
(Tj  can  be  determined  from  the  magnitude  of  (B-3),  since  the  ay  are  now  known. 
An  analogous  procedure  can  be  used  to  compute  /  o  and  ct.^. 


The  important  result  is  that  the  phase  angle  of  equation  (B-3)  is  proportional 
to  the  desired  "low  altitude"  Doppler  velocity.  The  j)hase  angle,  however,  will  be 
incorrectly  computed  if  the  amplitude  aj  is  negative.  In  order  that  aj  be  positive, 
it  can  be  shown  that  R  i{0)  must  be  greater  than  Wy  o  R-i^)/i>’o  o  or  about  0.6  Ro{Q). 
Owing  to  the  stochastic  nature  of  weather  echoes  and  the  limited  number  of  sam¬ 
ples  available  from  an  ASR  for  integration,  this  condition  may  not  be  met,  partic¬ 
ularly  if  the  autocorrelation  lags  for  the  two  beams  are  computed  on  alternate 
antenna  scans.  To  ensure  valid  determination  of  the  phase  angle,  we  solve  a 
modified  set  of  equations  where  the  high  beam  data  samples  are  scaled  to  have  the 
same  integrated  power  as  the  lf)w  beam.  The  resulting  estimate  for  the  "low- 
altitude"  mean  Doj'pler  frequency  /,  is  equation  (3)  in  Section  lI-B-2. 
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APPENDIX  C:  SIMULATION  OF  AIRPORT  SURVEILLANCE 
RADAR  SIGNALS  FROM  MICROBURSTS 


The  power  spectrum,  5,  measured  in  an  range-azimuth  cell  by  a  fan  beam  ASR 
can  be  expressed  in  terms  of  the  elevation  angle  resolved  field  of  velocity  spectra, 
5,  as: 

IT 

J 

fS(9,<f>,R,v)BTR{9)de 

=  ; -  (C-1) 

jBrA9)d9 

0 

where  Btr{9)  is  the  two-way  elevation  power  pattern  of  the  ASR  antenna.  Given 
either  measured  or  assumed  vertical  profiles  of  weather  velocity  spectra,  this  rela¬ 
tionship  can  be  used  to  synthesize  the  power  spectrum  that  would  be  measured  by 
an  ASR. 

As  an  example.  Figure  C-1  shows  power  spectra  that  would  be  measured  in  the 
model  "microburst"  described  in  Section  III-A.  Plots  in  the  left  column  are  for  the 
approaching  radial  velocity  cores  and  those  in  the  right  column  are  for  the  reced¬ 
ing  cores.  Center  ranges  of  3,  6,  9  and  12  km  are  assumed.  Recall  that  in  both 
cores,  the  magnitude  of  the  surface  radial  wind  component  is  taken  as  15  m/s 
(dashed  vertical  line),  decreasing  to  5  m/s  with  opposite  sign  1000  m  above  tne 
surface.  The  simulations  are  in  good  agreement  with  the  bimodal  spectrum  model 
assumed  in  this  report.  The  more  complex  structure  sometimes  observed  in  meas¬ 
ured  spectra  results  from  vertical  wind  profiles  more  complicated  than  assumed,  or 
from  statistical  error  in  the  spectrum  estimates  (see  below). 

The  stochastic  nature  of  radar  signals  scattered  from  precipitation  can  be 
simulated  using  a  Monte  Carlo  method  proposed  by  Zrnic  (1975)  and  used  by  Sir- 
mans  and  Bumgarner  (1975).  Equation  (C-1)  is  evaluated  on  a  discrete  grid  of 
equispaced  frequencies  spanning  the  Nyquist  interval.  A  single  realization  of  the 
discrete  Fourier  transform  (DFT)  of  a  signal  conforming  to  this  spectrum  shape  is 
then  simulated  by  multiplying  the  square  root  of  the  spectral  lines  by  randomly 
generated  complex  numbers.  The  amplitudes  of  these  numbers  are  Rayleigh  dis¬ 
tributed  and  the  phases  are  distributed  uniformly  between  0  and  27r.  An  inverse 
DFT  then  provides  synthetic  in-phase  and  quadrature  radar  signals  with  the 
appropriate  spectral  distribution.  Our  simulations  employed  64-point  Fourier 
transforms,  from  which  the  34  samples  required  for  emulation  of  the  signal  pro¬ 
cessing  operations  described  in  the  text  were  extracted. 

Simulation  of  spatial  smoothing  was  accomplished  by  generating  the  appropri¬ 
ate  number  of  independent  signal  realizations.  Adjacent  range  and  azimuth  gates 
in  real  radar  data  are  not  fully  independent  because  radar  pulse  and  beam  shape 
create  an  overlap  region  where  the  same  scatterers  contribute  to  the  echoes.  This 
effect  was  not  simulated  and  would  result  in  estimate  variance  slightly  larger  than 
our  calculations. 

Figures  C-2  are  simulated  power  spectral  estimates  for  the  underlying  spectra 
in  Figure  C-1.  Incoherent  averaging  of  three  independent  realizations  has  been 
performed  to  reproduce  the  range  averaging  applied  to  our  real  data.  Substantial 
distortion  of  the  underlying  shapes  occurs  owing  to  the  large  variance  for  power 


Power 


spectrum  estimates  with  three  degrees  of  freedom.  In  general  however,  the  esti¬ 
mate  fidelity  is  sufficient  to  show  the  overall  width  of  the  spectra  and  the  impor¬ 
tant  low-high  beam  amplitude  relationship  in  the  velocity  interval  associated  with 
near  surface  scattering. 

Signals  from  an  ASR-9’s  8/10  pulse  alternating  PRF  waveform  can  be  simu¬ 
lated  by  expanding  the  frequency  domain  for  spectrum  synthesis  to  a  value  equal 
to  the  reciprocal  of  the  largest  time  increment  that  is  a  submultiple  of  both  pulse 
repetition  intervals.  After  inverse  Fourier  transformation,  the  non-equispaced 
data  samples  are  extracted  from  the  resulting  oversampled  sequence.  Using  simu¬ 
lated  alternating  PRF  data.  Figure  C-3  repeats  the  calculations  of  velocity  esti¬ 
mate  bias  and  standard  deviation  shown  previously  in  Figure  III-l.  As  described 
by  Weber  (1987),  a  four-coefficient  shift-variant  interpolation  filter  has  been  used 
to  reconstruct  a  uniformly  sampled  data  sequence  prior  to  clutter  filtering  and/or 
velocity  estimation.  The  bias  of  velocity  estimates  would  be  unchanged  by  use  of 
the  alternating  PRF  waveform;  estimate  standard  deviation  increases  slightly  rela¬ 
tive  to  the  constant  PRF  calculation.  After  spatial  smoothing,  calculated  esti¬ 
mate  standard  deviation  varies  from  0..b  to  1.0  m/s  increasing  with  range  and  the 
magnitude  of  the  true  surface  velocity.  As  previously,  larger  standard  deviations 
apply  for  low  Doppler  weather  when  a  clutter  filter  is  used  (part  c). 

Note  that  the  interpolation  procedure  would  break  down  for  weather  signals 
that  exceed  the  Nyquist  velocity  associated  with  the  lower  PRF  (about  25  m/s). 
In  this  situation,  aliasing  to  different  parts  of  the  velocity  spectrum  would  occur 
at  the  two  PRFs,  causing  potentially  significant  spectrum  distortion.  To  prevent 
resulting  velocity  estimate  errors,  tests  should  be  performed  on  received  signal 
parameters  within  each  PRF  block  (for  example,  the  first  spectral  moment)  to 
detect  differential  aliasing.  Appropriate  dealiasing  procedures  couhl  then  be 
applied  to  the  signals  prior  to  velocity  estimation. 
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SPECTRAL  DIFFERENCING 


AUTOCORRELATION  BASED 
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(a)  Range  of  6  km  and  no  clutter  filtering  are  assumed. 


Figure  C-3.  Spectral  differencing  and  autocorrelation  based  velocity  estimate  bias 
and  standard  deviation  versus  "true"  outflow  velocity.  The  simulations  used  here 
reproduce  the  alternating  PRF  waveform  of  an  ASR-9  and  the  use  of  a  four- 
coefficient  interpolation  filter.  Plot  format  is  as  in  Figure  III-l. 
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